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1. Introduction 

There are a few exact solutions of the Einstein equations describing the gravitational 
field inside a stationary and axisymmetric rotating perfect fluid, the basic candidates to 
form a stellar model in General Relativity [1, 2]. Among them, only one is known to admit 
a closed surface of zero pressure, the key component to build a stellar model matching the 
interior (source) spacetime with a suitable asymptotically flat exterior. It is the Wahlquist 
metric, that describes a rigidly rotating perfect fluid with a 4-velocity expressible as linear 
combination of the Killing vectors, possesses the mass density-pressure equation of state 
(EOS) fj, + 3p = fJLo and has Petrov type D [3]. Nevertheless, it has been shown in several 
different ways that it can not correspond to an isolated object nor be matched with an 
asymptotically flat exterior [3-5] . Accordingly, General Relativity still lacks of any exact 
solution that can describe the interior of such stellar model. 

To find these global models, numerical methods and analytic approximations have 
therefore been the way to go. Some of us introduced a new analytic approximation scheme 
in [6, 7] focused in this kind of problem. It is a double approximation. The first one is 
post-Minkowskian with associated parameter A, which is related with the strength of the 
gravitational field and the second one is a slow rotation approximation with parameter 
fi, measuring the deformation of the matching surface due to the rotation of the fluid. 
We applied this scheme to find an approximate global solution to the gravitational field 
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of a fluid with simple barotropic EOS in some cases. We found solutions for constant 
density in [7], for a polytropic fluid (see [8]) and more recently for a lineal equation of 
state fi + (1 — n)p = [aq, [9-11] (hereinafter CGMR). When n = —2 this last case leads to 
the EOS of the Wahlquist solution. 

It is worth noting that despite Wahlquist's inherent interest as an exact solution, our 
metric is in some sense more general even fixing n = — 2. This comes from the fact that 
although Wahlquist is the most general Petrov type D solution for this EOS, symmetries 
and motion of the fluid [12], our metric can also be of type I. More interestingly, our 
solution, unlike Wahlquist, can be matched to an asymptotically flat exterior and thus 
can describe a compact isolated object. 

Another problem of Wahlquist's solution is its rotation parameter. Since Wahlquist 
is written in a co-rotating coordinate system, the explicit relation between the timelike 
Killing field and the rotation is not known. All that is known is that letting the Wahlquist's 
tq parameter go to zero in a limiting procedure (that involves a coordinate change that is 
singular when ro —> 0) we get the static spherically symmetric Whittaker solution [2, 3]. 
Despite this singular character of this limit, the relation it shows between Whittaker 
and Wahlquist seems to be a sound one since in the slow rotation formalism of [13], the 
Whittaker- like and Wahlquist-like slowly rotating metrics coincide [14, 15]. 

The question we try to solve in this paper is whether or not we can include an appro- 
priate approximation of the Wahlquist solution in our family of approximate solutions. 
There are two ways to answer. One of them is asking our n = —2 solution to be of Petrov 
type D, since a metric with its characteristics and this Petrov type must belong to the 
Wahlquist family [12]. The other one is finding a coordinate change to make them to 
coincide. 

Regarding the first one, we have already verified that our solution can take Petrov 
type D in [10, 11, 16]. Some of its free constants — the ones unrelated with gauge at first 
order — are then fixed and we found their values do not coincide with the ones they are 
forced to take when we matched our interior solution with an asymptotically flat exterior 
solution, as one expects. 

The coordinate change way involves writing our solution in a co-rotating frame, and 
then making a series expansion of Wahlquist's solution with /j,q as post-Minkowskian 
parameter — i. e. when this parameter equals zero the Wahlquist solution becomes 
Minkowski's metric, the same behaviour parameter A has in our solution — . This way is 
more meaningful since, in spite of any result we can get from our approximate metric alone, 
there is always the question of whether our solution really corresponds to a parametric 
expansion of an exact metric. Working this way we verify explicitly this correspondence. 

In Section 2 we give some notation and definitions used along the paper and we write 
the CGMR metric for n = — 2 and perform the rotation to write our metric in a co-rotating 
coordinate system. In Section 3 we give the Wahlquist metric, write it in spheroidal-like 
coordinates and then write the approximate post-Minkowskian Wahlquist metric. Finally, 
in Section 4 we compare both solutions and determine the value of our constants and the 
relation between ro and the rotation parameters. 



2. The approximate interior metric 



We work within the analytical approximation scheme developed in [6, 7] and [11]. It 
allows to build an approximate stationary and axisymmetric solution of the Einstein's 
equations for a source spacetime and an asymptotically flat vacuum region around it, 
although in this paper we put the focus on the interior, source spacetime. This section is 
devoted to a brief review of its main points in the general formulation. 

Let's choose t and (p to be the coordinates adapted to the Killing vectors £ and £, and 
r and 6 to the orthogonal two-dimensional surfaces [17, 18], then we can write the interior 
-and exterior- metric with the structure 



in the associated cobasis u> = dt, u r = dr, or = rdO, = rsinOdip. We require these 
coordinates to be spherical associated to harmonic Cartesian- like ones. 

The non convective 4- velocity can be written as it = ip(£ + Here, ip is a normal- 
ization factor and uj the angular velocity of the fluid. The explicit expressions for \x and 
p for the linear EOS with n = —2 are 



with ip-£ the value of ip on the p = surface. 

For n = —2 the metric CGMR obtained in [11] is the following (we have made the 
notation change: tq —> r s to avoid confusion with the ro parameter that appears in the 
usual expressions of the Wahlquist metric): 
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where, to simplify we have taken the gauge constants ao = 0, a 2 = 0, 6 2 = because they 
are not needed hereafter. The approximation parameters have the expressions 
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and 5 appears in the series expansion of ^s, i.e., = 1 + A5. Its worth noting ft gives 
a measure of the deformation of the source due to its rotation. This solution is apparently 
less interesting than the Wahlquist exact solution for the same kind of source because it is 
an approximation. Nevertheless, it is more general in a sense because it is a Petrov type 
I solution unless 



m 2 = ^ + 0(A,ft 2 ), j 3 = ^ + 0(A,ft 2 ), 
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in which case it becomes a Petrov type D solution [11]. Moreover, in our solution the 
parameter uj has the clear meaning of angular velocity of the fluid, 



(11) 



and its vanishing leads to a static solution (i. e. jttp = 0) while there is no equivalent 
parameter in Wahlquist. Another interesting property of our solution is that, in general, 
it can be matched to an exterior asymptotically flat solution which is not the case for 
Wahlquist metric. 

If we want to compare this approximate solution with Wahlquist 's metric we have to 
start finding their expressions in the same coordinates. The first problem is that Wahlquist 
metric is written in a co-rotating coordinate system and CGMR is not, so first we must 
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choose between the two kinds of coordinates. Changing CGMR to a co-rotating system 
is straightforward doing 
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and then in the co-rotating system the metric components are: 
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and the other components remain unchanged. Let us remark that the 7^ component is 
now of order A 1 / 2 instead of the order A 3//2 it was in the original coordinates (see [11] for 
some comments). 



3. The Wahlquist metric 

The next steps in the comparison are, using the singularity free Wahlquist metric, first 
expand it in the appropriate approximation parameters and then make coordinate changes 
to reduce it to a particular case of CGMR. 

We start writing the 4- velocity of the fluid as 

u = f- 1/2 d t (g al su a ^ = -1) , (15) 
the energy density and pressure are given by 

M = ^o(36 2 / - 1) 

(16) 

and the singularity free Wahlquist metric reads [2] 
ds 2 = -f(dt + Adip) 2 + 
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Here fiQ, b, r$ are arbitrary constants and 770 and c are adjusted so that the solution behaves 
properly on the axis. The symmetry axis is located at 77 = 770 where 



h2(vo) = 0, 

and to satisfy the regularity condition c must be 



dh 2 
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Therefore 770 and c become functions of the constants /jo and 6. 

The constants b and /Jo are the values of the normalization factor f^ 1 ^ 2 and the en- 
ergy density on the matching surface of zero pressure. The ro parameter allows to get 
Wahlquist's static limit — Whittaker's metric [19] — making the change 



R 

{£, v} -> {R,x} ■ {£ = —,v = cos x} 



(24) 



and letting ro go to zero [3]. This suggests a connection between ro and Wahlquist's fluid 
rotation speed, although it is by no means necessarily true. Actually, to the best of our 
knowledge, there is no proved explicit expression for such rotation speed. The parameters 
{ro, Ho} will be the natural choice for us to expand Wahlquist metric, but first we must 
find the change to spherical-like coordinates. 



3.1. Wahlquist written in spherical-like coordinates 



Our approximate metric (3) to (6), (13) and (14) is written in "standard" spherical 
coordinates — in the sense that in when A = they are the usual ones — so we need to 
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find a consistent way to write the Wahlquist metric in a set of coordinates as close to ours 
as possible to begin with. 

In this regard we note first that if we put hq = in the Wahlquist metric (17) we 
obtain Minkowski in oblate spheroidal coordinates {£, r/}, whose coordinate lines are oblate 
confocal ellipses and confocal orthogonal hyperbolas. From these coordinates is easy to 
go to Kepler coordinates {R,x} changing £ = i?/ro,ry = cos^j where R represents the 
semi-minor axis of the ellipses, x the Kepler eccentric polar angle and ro the focal length. 
Finally we get standard spherical coordinates {r, 6} by changing 



R 2 + Tq sin x 



r sin 9, R cos x = r cos 



(25) 



Moreover, the limiting procedure (24) from Wahlquist to its static limit (Whittaker) 
has a similar form, in this case leading to Kepler-like coordinates. 

These considerations suggest to look for a change of coordinates in the Wahlquist 
metric (prior to any limit) so that the new coordinates "directly represent" spheroidal- 
like coordinates. Fortunately this approach arises when we plot the graphics of and 
^2(77) (Fig. 1.) We can see that these curves describe locally a hyperbolic cosine and a 
squared sine, respectively. Taking all this into account we write 
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where we introduce R±, R2 just to simplify calculations later. Let us write now the two 
dimensional metric spanned by {£, 77} 

dY? = Ad£ 2 + Bdrf 
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FIG. 1. Behaviour of the and ^2(7?) functions for certain values of k and b. 
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taking hi, h 2 as functions of Ri and R 2 
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Now let us do another coordinate change to a kind of spherical coordinates {r, 6} using 
the previous relations Eqs. (25) and (26). Then, Ri and R 2 are the following functions of 
r and 9 
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Only the terms mn + 77722 and mn — 777-22 depend on both {/jo, 7"o}; the remaining terms 
depend on ro alone. 
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Now we are going to write the full metric in terms of {r, 9}; notice that the inversion 
can only be approximately done (in a series of no). First, we determine r]o up to order 
/Xq , and then c to the same order. This last series depend on b 2 and we want to compare 
with our approximation, where the pressure is of a higher order in the post-Minkowskian 
parameter (in this case in fio) than the density. To determine how b depends on no we 
recall that for no = the Wahlquist metric becomes Minkowski's metric written in oblate 
spheroidal coordinates. From Eq. (16) this means that 
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so if we want to have p = "P\n\ + C[Mo]> the constant must be b 2 = 1 + 0[mo]- Accordingly, 
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where a\ and o-i are two new constants. Using this hypothesis on b 2 we obtain for the 
constants t]q and c up to 0(/ig) 
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Next, we invert the change of coordinates Eq. (26), what gives 
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And finally, by doing the coordinate change {R\,R2} —> {r, 9} we obtain the metric 
coefficients up to C(Mo) i n the spherical-like coordinates desired 
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4. Comparing the approximate Wahlquist solution with the CGMR solution in 

co-rotating coordinates 

Now we face the problem of identification of the parameters and to perform the final 
adjustments of coordinates needed to make every term in Wahlquist and CGMR equal. 

To get an idea of the problems arising, we analyze first the static limit. By using the 
relation between mass parameters A = g/xo 7 " 2 an d making ro = in Eq. (42) we obtain 
the expression for the j rr coefficient of the static metric 
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and upon comparison with the corresponding static limit of CGMR 

cOMR, n \ , r 2 X 2r 4 A 2 2m r 2 X 2 24r 2 S\ 2 
7r c r GMR (f} = 0) = 1 + m A =- - — -j + — °— 9 — — + 0(A 3 ) (49) 



we can see that there are discrepancies among r 4 terms. To some extent this was to be 
expected since CGMR was written in coordinates associated to harmonic ones and no such 
a condition has been imposed on Wahlquist. In this particular case, the two metrics can 
be rendered exactly equal with a change of the form r — > r[l -h A 2 /(r)] . 



4.1. Adjusting parameters 



We go back now to the non-static case. If we compare the lowest order term in gt^ and 
gu of both solutions we can see that the constants A and f2 of the CGMR solution must 
be related with the and tq constants of Wahlquist as follow 
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where k gives the proportionality between the rotation parameters of both solutions, ro 
and f2, and will be determined later on. If we perform this identification we get to a new 
difficulty because the Wahlquist solution has A-free terms with f2 dependence. These terms 
appear associated with powers of /io^o [ or ^ 2 ^ 2 using Eq. (50)] . This is not possible in our 
self-gravitating solution building scheme. This issue can be solved using the remaining 
freedom in time scale and {r, 9} coordinates. The changes we can do are 
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with Fi constants depending on the parameters and Gj, Hi undetermined functions. Im- 
posing vanishing of these unwanted terms, we get the time scale change 
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and the {r, 9} changes 
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Note that the symmetry axis for the old coordinates is located at 6 = 0, ir and due to the 
presence of the sinO it remains at G = 0,7r. We will maintain this condition for all the 
coordinate changes of the 9 coordinate. 

Now, we introduce these changes in our last expression of Wahlquist metric obtaining 
up to C(/Uq) 
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After dealing with hq and ro, we have to find expressions for b and /c. Recalling Eq. (37), 
we wrote b 2 as a series in /j,q and o"i, <T2. However, according to the approximation, b 2 
should also contain terms in ixqTq. We get this behaviour if we write 
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o- 2 -> (o- 2 r 2 + rgi/ 2 )r 2 
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with the distance r s added for dimensional reasons (also appears in the change of constants 
(50) and in the CGMR solution) which represents the radius of the non-rotating matching 
surface. In terms of A and Q, this leads to 
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with a% and fii constants to be adjusted. 

Now we can write the approximate Wahlquist metric in terms of our parameters A and 
Q using (50). Comparing the lower terms in A for gt v of the CGMR co-rotating solution 
and the approximate Wahlquist solution just built, we can determine the proportionality 
constant k to be a series in our rotation parameter SI 

K = l-^ + 0(n 3 ) (65) 



4.2. Adjusting terms 

Once the relations between the approximation parameters of both metrics are deter- 
mined we can obtain the expression of the approximate Wahlquist metric written in the 
same parameters we have used for the CGMR co-rotating solution. To make both so- 
lutions coincide the remaining freedom is a change in the {r, 8} coordinates of the type 
displayed in Eq. (52). If we make this change in the Wahlquist metric 



r < 1 + AST 3o"i sin 



lr 2 
2^! 



cos 



y r 
5 r 



2 <71 + 



2r 4 
7ri 



+ sinflcosfl 



70 
lr 2 



1 A ,r 4 /13 o 
n ~i +33 cos 2 
rf V 3 



+ o(A 3 , n 4 ), 



2, 2+3 ^-2T0 A ^) + °( A3 '^ 



(66) 
(67) 



we get that, for the two metrics to be exactly equal up to C(A 2 , Q 3 ) the free parameters 
of the CGMR co-rotating solution must be 



m = 0(A 2 , 

j x = ®n 2 s + o(\, o 4 ) 

5 



m 2 = |(1 + 2AS) + 0(A 2 , ft 2 ) 
5 



33 



36 
175 



+ c(A, n 2 



(68) 
(69) 



and the free parameters of the approximate Wahlquist solution must be 



ai = S, fj 2 = 0, isi 



»2 
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(70) 



This gives b 2 as 



[i + n 2 )(i + 2\s) + o(\ 2 , n A ). 



(71) 



This result agrees with the meaning of S and b. The term (1 + SI 2 ) comes from the change 
of the normalization factor over the transformation of the temporal coordinate Eq. (53) 
we have done . 
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The final expressions for the metric components in the orthonormal basis are, 0(A 3 , ft 4 ) 



Ire 



lee 



l w 



lt v 



Itt 



1-A 

+ ft 

31 

315 



ft 2 P 2 + ^A 



7r 2 



+ 65 



5 r 2 
3 r 2 



125 

„2 



lr 2 
7 r 2 



1-A- 



+ ft z 



1 - A- 



+ ft 2 



-ft 2 p 2 



21 



-A 



189 



5 



12 r 2 

5 r 2 



25 r 2 



+ 



1 



/181 r 2 
6 r 2 ' 5 V 3 r 2 



1265J P 2 



6 



ft 2 p 2 



105' 



-A 



1895- 12- 



17 r 2 
6 r 2 



110 r 



3 r 2 



~ - 1265 ] Pi 



-X^Sl- \ P\ + 



„2 



-P/ - 3ft 2 f 35 + 



2r' 
5r! 



35 



-1 - A- 



+ 



A 2 ^ 
5 r 2 



-ft 2 



+ ft 2 



1 + -P 2 
5 



8r 2 

3 r 2 



+ ( 125 



34 r 2 
21 r 2 



(72) 
(73) 



(74) 



(75) 
(76) 



(77) 



To verify the whole procedure we can compare now with the conditions necessary for 
our n = —2 approximate metric to be of type Petrov D [11, 16, 20], i.e., Eq. (10). They are 
compatible with the values of the constants m 2 and j'3 we have just found in (68) and (69), 
as wished. Also, when matched with an asymptotically flat exterior, m 2 , j'3 and the rest of 
the metric free constants can only have the expressions we found in [11]. Since the n = —2 
fluid for a type D interior does not satisfy the matched expressions, we concluded then 
that it can not be the source of such exterior. Nevertheless, it is worth noting here that 
CGMR contains a n = —2 sub-case that lacks this problem and can indeed be matched 
that way. It has then all the characteristics of Wahlquist's fluid but Petrov type I instead 
of D. 

Note, finally, that the Cartesian coordinates associated to the spherical-like coordinates 
used above are not harmonic. Nevertheless, since Eqs. (72) to (77) correspond as well to 
the stopped n = —2 CGMR with particular values of the free constants, undoing the 
change (12) they become harmonic again. 
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5. Remarks 

In this work we have taken the singularity free Wahlquist metric and managed to trans- 
form it into the form the CGMR metric takes when written in a co-rotating coordinate 
system. 

We have identified Wahlquist's parameters corresponding to A and f2 of [6]. Doing this, 
we have found an expansion of the parameter ro of Wahlquist's metric in terms of our £1. 
Accordingly, now we have an approximate expression of ro in terms of physical quantities 

• s -^ n ( i -w) +o(n4) -s w ( i -^) +o (S) (8 ' ^' (78) 

To the best of our knowledge its qualitative relation with the rotation was only guessed 
through the singular limiting procedure that takes Wahlquist and leads to Whittaker's 
metric but no parametrization of it in terms of well defined quantities had been given. 

Last, notice that the usual interpretation of uj = ^ as rotation speed of the fluid 
as seen from the infinite lacks sense if we deal with a metric that is not matched with 
an asymptotically flat exterior. Besides, the definition of stationarity and axisymmetry 
allows a change of coordinates {t = t', <p = ip 1 + at'} that can modify the value of uj to 
uj' = uj — a or make it zero (the case of co-rotating frames). Nevertheless, when dealing 
with a family of metrics explicitly dependent on uj, its value can be important. In the case 
of, e.g., CGMR, we see that written in co-rotating coordinates Ut'/u^i = but uj is part 
of the metric functions and actually, uj — > still leads to a static metric. In this sense, the 
characterization of ro (78) is meaningful. 
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